GUTPA/03/04/01 



ON 
(N 



Supersymmetric Grand Unification 
with a Fourth Generation? 



cn ' J. E. Dubicki f and C. D. Froggatt t 

o ■ 
o 

C\| \ Department of Physics and Astronomy 

' University of Glasgow, Glasgow G12 8QQ, Scotland 

Abstract 



The possibility of incorporating a chiral fourth generation into a 
SUSY-GUT is investigated. Precision fits to electroweak observables 
' require us to introduce light supersymmetric particles, with masses 

less than Mz- These particles might also provide decay channels for 
. the fourth generation quarks of mass ~ 100 GeV. We also require tan j3 

to lie in the range 1.50 < tan/3 < 1.75 and obtain an upper limit on 
the lightest Higgs boson mass in the MSSM4 of 152 GeV. 

o 

*§h' 1 Introduction 

The Standard Model (SM) is, of course, phenomenologically very successful, 
while its supersymmetric extension stabilizes the gauge hierarchy problem 
and allows a grand unification of the SM interactions. However, there is 
no explanation of why there should be just three generations of quarks and 
leptons or their hierarchy of masses. In this paper we investigate the possibil- 
ity of consistently incorporating a fourth generation into a supersymmetric 
grand unified theory. 

We shall assume a structure akin to that of the minimal supersymmetric 
standard model (MSSM3), adding a complete chiral fourth generation and 
its associated supersymmetric partners (the so-called MSSM4). We note that 
at one-loop level the validity of gauge coupling unification is independent of 
the number of generations, and so this requirement does not discriminate 
between three and four generation models. However, the requirement that 
all the Yukawa couplings remain perturbative at energies up to the GUT 
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scale places restrictive upper bounds on the fourth generation T, B, E and 
N masses and also constrains the supersymmetric parameter tan j3 = ^ 
{vjj and Vd being the vacuum expectation values of the two Higgs doublets 
Hjj and Ho that are present in the MSSM). In section |21 we investigate 
the masses and decay channels of the fourth generation particles that are 
consistent both with perturbative gauge coupling unification and the latest 
experimental bounds from direct searches. In section |3] we comment on the 
issues regarding fits to precision data in the MSSM4 before discussing the 
lightest Higgs boson mass in section 0] 

2 Experimental Constraints and a Fourth Gen- 
eration 

All fourth generation models must adhere to certain experimental constraints, 
the first of which stems from precise measurements of the decay charac- 
teristics of the Z-boson performed at LEP. This has set a lower bound of 
Mp > on an y non-SM particles that couple to the Z-boson. Ignoring the 
unnatural hierarchy emerging within the neutrino sector, we assume a Dirac 
mass ~ i^F) f° r the heavy neutral lepton. 

We begin our discussion by considering the leptonic sector where we as- 
sume Me > M^. Under the assumption that the mixing between the fourth 
generation leptonic sector and the first three generations is negligible, the 
decay E — ► NW* will be dominant. Current experimental limits searching 
for E — > NW* from e + e~ — > E + E~ production have been performed by the 
OPAL and L3 collaborations up to the kinematic limit Me ~ 100 GeV pQ. 
However, it turns out that to be consistent with perturbative unification we 
require Me ~ so that the fourth generation charged lepton Yukawa cou- 
pling Ye doesn't become non-perturbative below the GUT scale. Therefore, 
in order to evade experimental bounds, the mass difference AMl must be less 
than ~ 5 GeV, which is in the region where the trigger efficiencies are signif- 
icantly lowered and events are dominated by the two photon background [T] . 
We also note that with AM^ the order of a few GeV the decay lifetime 
t (E —+ NW*) is too short for the heavy E lepton to leave a charged track in 
the electromagnetic calorimeter. Regarding the heavy neutrino, OPAL and 
L3 have set the bound M N > 70 - 90 GeV based on the search for iV — > IW* 
(I = e, fx or r) provided the mixing matrix elements satisfy Vjv eiMiT > 10 -6 p. 
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However, if we assume that this mixing angle is negligible {Vjq e ^ T < 10~ 6 ), 
then the neutrino is stable enough to leave the detector. In this case, the 
only relevant bound comes from the Z° decay width noted earlier. Based on 
the above discussion, in the rest of this paper we shall take: 

M £ ~50GeV ; M N ~ 50 GeV (1) 

as representative masses of the fourth generation leptonic particles. 

Direct searches for the fourth generation quarks is an ongoing process at 
the Fermilab Tevatron. Here we focus on the experimental restrictions for 
M t > Mt > Mb, so that the charged current (CC) decays B — > tW~ and 
B — > TW~ are kinematically forbidden. The leading CC decay mode will 
then be B — > cW~ which is doubly Cabibbo suppressed by the mixing matrix 
factor V cB . In this situation, loop induced flavour changing neutral current 
(FCNC) decays can dominate provided |2]: 

1^1 < O (10- 2 - 10- 3 ) (2) 

I V tB\ 

Several experiments have searched explicitly for B quarks decaying via FC- 
NCs. The .D0 collaboration [3. has excluded the range < M B < 

Mzo + Mb by a null search for both B — > try and B — > bg. For masses 
M t , M T > M B > M z o + M b , the decay B — > bZ° is expected to dominate, 
except for B -> bh°, if M B > M h o + M b . The CDF collaboration @j has per- 
formed a general search for long-lived particles that decay into a Z° gauge 
boson. This will encompass the FCNC decay of a fourth generation B quark 
decaying via B — > bZ°, if the mixing matrix factor VtB — Vrb is small enough 
to result in a long lifetime. By looking for Z° — ► e + e~ with a displaced vertex 
they are able to exclude a B quark mass up to M B = 148 GeV for cr B = 
1 cm, where t b is the proper decay time of the B quark, and a branching 
ratio of Br (B — > bZ°) = 100 %. However, this limit diminishes to M B ~ 96 
GeV if ct b > 22 cm. To date, this remains the only lower mass bound on 
quasi-stable B quarks as emphasized by Frampton et al. [3]. 

As regards the T quark, for Mt ^100 GeV, there are two competing decay 
modes; T — > BW and T —>■ bW. The BW decay will certainly dominate over 
the bW decay when the former is real, since the bW channel is suppressed by 
the mixing matrix factor Vrb- However, when Mt~M b < Mw, the two body 
bW decay could be competitive with the three body BW* decay. For large 
enough Vrb the T — > bW decay will be dominant. However, for Mr ^ M t , 
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such a dominance would lead to a large excess of bW events relative to those 
already present from tt production and decay, and so this implies that 
must be small enough for the T — ► BW* decay to dominate. The detection of 
the T quark would then depend on the decay properties of the B quark and 
the mass difference Mt — Mb- We refer to [B] and [7j for a discussion of this 
scenario, in particular Gunion et al. jn] finds the only way that TT events 
can evade being included in the CDF and .D0 data sample is if Mt — Mb is 
sufficiently small so that the W* in T — > BW* is virtual and the jets and 
leptons from the two W*'s are soft. A further analysis on updated data is 
needed that takes into account the quasi-stable nature of the B quark. 

To obtain the allowed masses of the fourth generation quarks that are 
consistent with perturbative gauge coupling unification we perform a renor- 
malization group study of the MSSM4. Specifically, this enables us to place 
upper limits on the masses of the T and B quarks by ensuring their Yukawa 
couplings Y T and Y B run perturbatively to the GUT scale M GUT : 

Yt,b O) < 4?r f o r M zo < fJ. < M GUT (3) 

where Mqut is defined to be the scale where ai(/i) = 0:2 (/-0- In our analysis 
we have neglected all the Yukawa couplings from the first three generations, 
except that of the third generation t quark whose mass we take to be M t = 
175 GeV. As is typical with four generation models, we also require small 
values of tan/3 so as to avoid Y B (M z o) > O (\A47r). Further details of the 
procedure for running the renormalization group equations in the MSSM4 
can be found in jHj . In Figure ^ we have plotted the maximum T quark mass 
(M™ aa; ) versus Mb, given that perturbative unification must occur. We plot 
for tan/5 = 1.60, 1.65 and 1.70, and we fix M t = 175 GeV. We have taken 
Mb > 96 GeV, which corresponds to the absolute experimental lower bound 
on a quasi-stable B quark as discussed earlier. We can see from Figure [T] 
that higher tan/5 values allow for a larger M™ ax , though at the expense of 
restricting M B , and as such tends to favour the hierarchy: 

M T > M B (4) 

This mass hierarchy is necessary if we consider the T and B quarks to have 
the standard model decay channels T — > BW* and B — ► bZ°, as discussed 
earlier in this section . As we decrease tan/3 then the hierarchy in Eq. 0] 

1 If Mt < Mb then we would expect T — * bW to be the dominant decay channel of the 
T quark, which is excluded experimentally based on searches for the third generation top 
quark. 
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Figure 1: Plot of M^ ax versus M B for tan/? = 1.60, 1.65 and 1.70, as obtained 
from the requirement of perturbative gauge coupling unification. The third gen- 
eration t quark mass is taken as Mf = 175 GeV. The fourth generation leptonic 
masses are taken as Me = Afjv = 50 GeV. 

125 | 1 1 1 1 1 1 




96 98 100 102 104 106 108 

M B (GeV) 

Figure 2: Same as Figured but with M t = 170 GeV. 
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becomes harder to maintain. Overall, we find that tan [3 can take on the 
values 1.50 < tan/9 < 1.75 without the Yukawa couplings becoming non- 
perturbative below the GUT scale. In Figure El we also plot M™ ax versus 
Mb, but this time using M t = 170 GeV. We can see that this value of M t 
increases M™ ax by a few GeV for a given Mb and tan /3. For M t = 180 GeV 
the allowed ranges of Mt and Mb are very constrained, in fact we only find 
solutions for tan/3 > 1.66 and M T , M B < 100 GeV. 

Looking at Figures [T] and 121 it is clear that we could arrange Mt and Mb 
so that the standard model decay channels: 

T -> BW* ; B — > bZ° (5) 

are kinematically accessible, and dominate on the basis that the mixing angles 
Vbc and Vrb are suppressed. An updated analysis based on RUNII data 
at Fermilab is essential to exclude this scenario. However, we also note 
that there are other decay channels for the T and B quarks that become 
possible in the MSSM4. For instance, we can consider the possibility of light 
super symmetric particles providing decay channels for the T and B quarks. 
In this situation one can constrain the masses of the light (i.e. < M z o) 
neutralino (xl) and chargino (xt) pair which, as we will see in section El is 
already required by fits to precision data, in order to allow the following two 
body decays: 

T^B x xt ; B^B lX t (6) 

where Bi is the lightest mass eigenvalue in the fourth generation sbottom 
sector. It has already been shown by Gunion et. al. [6j that B\ is typically 
the lightest squark in the MSSM4. Ensuring both of the decay channels 
in Eq. |U] are kinematically accessible, combined with the constraints from 
perturbative unification and precision data fits, places severe restrictions on 
the allowed spectrum. We discuss this in more detail in section 01 

3 Precision Measurements and a Fourth Gen- 
eration 

It is difficult to provide bounds from precision data without a fully consistent 
study taking into account exact particle masses, any light supersymmetric 
particle spectra present and mixings between different flavours. However, 
Maltoni et. al. [9 :j has pointed out that a highly degenerate neutralino and 
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chargino pair can provide the necessary contributions to the precision param- 
eters that are needed to cancel that of the fourth generation, whilst at the 
same time being consistent with LEP bounds [TUI . Specifically they require: 



M w+ < AM X < 3 GeV and M x ~ 60 GeV (7) 

where we define the notation AM X = M-± — M^o and M x = M-± > M^o. 
Looking at their results, we see that the magnitude of the contribution to 
the fitted parameters from this sector is highly dependent on M x . Deviations 
from My, larger than + 30 GeV or - 5 GeV are ruled out at the 2a level. We 
now discuss the influence of these precision data fits on the scenario where 
the T and B quarks decay into supersymmetric particles. 

If allowed, the B — ► B%Xi decay in Eq. [Hjwill certainly dominate over the 
one-loop FCNC decay B — > bZ° and two-generation decays B — ► cW~ that 
traditional searches have looked for. Typical (pole) masses required for this 
scenario to be viable are: 

Mb (— Mt) > 100 GeV 

M Bi < 50 GeV (8) 
55 GeV < M^ < M B - M Bi 

We assume that 9 B ~ 1.17 rad, where 6 B is the mixing angle that relates the 
Bl and Br weak eigenstates to the mass eigenstates B\ and B2 through the 

B 2 J I -sm9 B cos9 B J I B R 



(9) 



This value of 8 B is chosen so that B\ decouples from the Z°-boson, thereby 

f- < M z< 



maintaining consistency with the measured total Z° width if M B < It 



turns out that 9 B ~ 1.17 rad is also required in order to evade the experi- 
mental direct searches for the light Bi squark. We discuss this in more detail 
later in this section. 

The decay rates for T — > Bix{ and B — > B\x\ depend crucially on the 
amount of Bl — Br mixing. It is especially important to consider the effect 
of fixing 9 B ~ 1.17 rad on the T — > B\xt decay. The Lagrangian for the 
B\-T-Xi interaction is given by [TT] : 

^TB^xt = gT P R (-U n cos 9 B + U 12 Y B sin 9 B )xfBi 

+gTP L {Y T V l2 cos 9 B ) xfBi + h.c. (10) 
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where P LtR = \ (1 =F 75) and Uy and Vij = 1,2) are the 2x2 unitary 
matrices diagonalizing the charged gaugino-higgsino matrix 2 : 



M 2 y/2 sin (3M W ( M~± 



~ 1 >/2cos/?A% ^ ; \ M-± ' 

We take fin and M 2 to be in the parameter regions where we expect a light 
degenerate chargino-neutralino pair, required in the MSSM4 by fits to preci- 
sion data p. If we take M 2 <C \hh\ then the Xi an d Xi 1 are both gaugino-like, 
with M~± ~ M 2 . In this case the mass of the lightest neutralino is given by 
M^o ~ min (M!,M 2 ), where Mi is the U(l) gaugino soft mass. Therefore, 
small mass splittings only occur if M 2 < M\. The mixing matrix elements in 
the chargino sector for the gaugino-like case are given by: 

V u ~ Un ~ 1 V 12 ~ U 12 ~ (12) 

On the other hand, if <C M 2 then %5 an d Xi are both higgsino-like 
with degenerate masses M^± ~ M^o ~ In this case, the mixing matrix 
elements in the chargino sector obey: 

Vu ~ U\\ ~ V12 ~ sg'n (//#) t/12 ~ 1 (13) 

We can see from Eq. II 01 that the decay rate for T — > Bix[ will be suppressed 
if the light chargino is gaugino-like, since V12 ~ ^12 ~ and cos#g is 
fixed at 0.39. Therefore, in order to ensure that T — > Bixt dominates over 
T — > bW + , without appealing to a suppression of the CKM mixing matrix 
element Vrb, we would expect the light chargino xt t° De higgsino-like. 

As with most phenomenological studies of supersymmetric models with 
_R-parity, we must ensure that the lightest supersymmetric particle (LSP) is 
neutral for cosmological reasons. This is usually taken to be the neutralino 
Xi in the MSSM3. The masses in Eq. |H] contradict this requirement. One 
possible solution presents itself, however, if we assume that the LSP is in 
fact the fourth generation sneutrino N± which becomes stable due to im- 
parity 3 . Since Mfi < we must arrange the mixing angle 9$ in the 



2 fiH is the bilinear term that couples the two Higgs doublets in the superpotential, 
whilst M2 is the SU (2) gaugino soft mass that appears in the supersymmetric breaking 
Lagrangian. 

3 Previous studies of the MSSM4 have shown that for some of the parameter space it is 
natural to assume that N\ is the LSP 0- 
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fourth generation sneutrino mass matrix such that N\ decouples from the Z°- 
boson. A right-handed LSP Ni does indeed not couple to the Z°-boson. The 
light fourth generation bottom squark will then decay via the semi-leptonic 
channel B\ — ► clN*, where I = e,fj,,r. Such a decay involves the factor 
IVbcVmI, leading to a long lifetime. Experiments looking for hadronizing 
sbottom squarks currently only exclude the range 5 GeV < M§ < 38 GeV, 
if the mixing angle 8^ ~ 1.17 rad such that B\ decouples from the Z°- 
boson ^2] • Searches for long lived charged particles in pp collisions at CDF 
have only excluded particles in the mass range 50 GeV - 270 GeV |13j . 

4 The Higgs Sector of the MSSM4 

At tree-level the lightest Higgs boson mass M^o in the MSSM4 is bounded 
from above by the relation: 

M h o < M\cos2f3\ < M z o 
M = min(M z o,M A o) (14) 

This has already been ruled out from experimental data taken at LEP [14 . 
It is therefore of vital importance to check that the radiative corrections in 
the MSSM4 can provide a large enough contribution to raise the mass of the 
lightest Higgs boson to be above its experimental lower bound. 

To take the radiative corrections into account in the MSSM4, we use the 
one-loop effective potential, given by: 

V eff = V + V 1 + ... (15) 

where Vq is the tree-level potential, V\ contains the one-loop contributions, 
and the ellipses represent the higher-loop corrections which we shall ignore. 
The one-loop radiative correction to the effective potential in the MSSM4 is 
given by: 

* = E sbf- 1 )"** 2 * + ( hl § - 1) < 16 » 

k \ ' s 

where the sum is taken over all particles in the loop; Ck = 6(2) for coloured 
(uncoloured) fermions; </& are the spins and are the field dependent masses 
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Figure 3: Plot of the lightest Higgs boson mass M h o versus Ms- The upper and 
lower solid curves correspond to My = Mb = 100 GeV ; Me = M^ = 50 GeV 
and tan/? = 1.65 for the maximal-mixing and no-mixing cases, respectively. The 
dotted curves correspond to the uncertainty in M h o from allowing Mt, Mb and 
tan (3 to take on any value as long as perturbative unification occurs. 

of the particles in the loops at the renormalization scale Q. We choose to 
minimize the potential at the scale Q = max(M t , Mt)- 

In this analysis we will consider contributions to V\ that arise from the t, 
T, B, E and N fermions and their corresponding superpartners. Since we 
are at low tan (5 we can ignore contributions to V\ from the third generation 
(s)bottom (s)quark. To calculate the lightest Higgs boson mass we make the 
approximations: 

M h M A0 ,M^ = O(M s ) ; l = E,N ; q = t,T,B (17) 

where Ms represents the scale of the soft supersymmetric breaking terms. 
Further details of the procedure can be found in 

In Figure El we plot the mass of the lightest Higgs boson mass M^o versus 
Ms for M t = 175 GeV, and values of Mt, Mb and tan/3 that results in 
perturbative unification. The bold curves correspond to taking Mt = Mb = 
100 GeV, M E = M n = 50 GeV and tan/3 = 1.65. We have also determined 
the allowed (M T , M B , tan/3) parameter space that results in perturbative 
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unification, and from these sets of values we calculate the lightest Higgs 
boson mass. For a given Ms, we retain the maximum and minimum values 
of Myfi returned, which are represented by the dotted curves. We have also 
plotted for the maximal- mixing and minimal-mixing scenarios, as discussed 
in detail by Espinosa 15J. Overall, we can see that the maximum Higgs 
boson mass M^ per consistent with perturbative unification has a value of 
Mapper ^1^2 GeV. This is safely above the LEP lower bound. 

5 Conclusions 

We have seen that is it possible to incorporate a fourth generation into a 
super symmetric GUT model, requiring the existence of a light, degenerate 
neutralino and chargino pair in order to provide the necessary cancellations in 
precision data fits. The fourth generation masses are then tightly constrained 
with typical values of — Mb — 100 GeV and M E ~ M N ~ 50 GeV. To 
retain perturbative consistency to the unification scale we constrain tan f3 to 
lie in the interval 1.50 < tan/3 < 1.75. In order to provide decay channels 
to the fourth generation quarks, it might be that the LSP is the sneutrino 
Ni with a mass M^ < 40 GeV. We would also expect a light B\ with 
M Bi ~ 40-50 GeV along with a degenerate neutralino and chargino pair 
M^ ~ 55-60 GeV. Such a light B\ squark is chosen to decouple from the 
Z°-boson, but would be copiously produced at the Fermilab Tevatron and 
should be searched for. The upper limit on the lightest Higgs boson mass 
is M^Q Per ~ 152 GeV. The supersymmetric spectrum needed to satisfy all 
these constraints cannot be obtained from MSUGRA scenarios with universal 
parameters at the unification scale. We should also remark that a large degree 
of fine-tuning of parameters is involved in such a model. 
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